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Two-dimensional  Strain  Cycling  in  Plasticity 

by 

P.  M.  Naghdi  and  0.  J.  Nikkei,  Jr. 


Department  of  Mechanical  Engineering 
University  of  California,  Berkeley,  CA  94720 


Abstract  Detailed  calculations  are  presented  for  strain  cycling  in  a 
homogeneous  deformation  that  can  be  sustained  by  a  biaxial  state  of  stress 
in  thin-walled  specimens  of  OFHC  copper.  These  calculations  are  made  with 
a  set  of  relatively  simple  constitutive  equations  within  the  framework  of 
the  strain-space  formulation  of  plasticity.  The  predicted  theoretical 
calculations,  carried  out  in  the  context  of  small  deformation,  are  in  good 
agreement  with  corresponding  available  experimental  results  for  saturation 
hardening  and  erasure  of  memory  in  two-dimensional  strain  cycling.  Also, 
with  the  use  of  the  calculated  results,  a  scalar  quantity  that  characterizes 
strain-hardening  is  plotted  as  a  function  of  plastic  strains.  Such  plots 
are  likely  to  be  useful  for  computational  purposes. 


1 .  Introduction 

In  an  interesting  series  of  experiments  pertaining  to  two-dimensional  strain 
cycling  in  plasticity  that  can  be  sustained  by  a  biaxial  state  of  stress  resulting 
from  combined  tension-compression  and  torsion  in  thin-walled  specimens  of  OFHC 
copper*,  Lamba  and  Sidebottom  (1978a, b)  observed  the  following  three  phenomena: 

a)  The  occurence  of  saturation  hardening  after  loading  from  an 
undeformed  state  and  cycling  along  a  strain  path  which  is 
essentially  an  ellipse  in  the  normal  strain-shear  strain  plane; 

b)  the  erasure  of  memory  after  the  material  has  reached  a  state  of 
saturation;  and 

c)  the  nature  of  the  post-saturation  stress  response  for  cycling 
in  a  relatively  "complex"  nonproportional  strain  path. 

With  reference  to  the  above  two-dimensional  strain  cycling  experiments,  our 
main  objective  is  to  examine  the  predictive  capability  of  a  (rate-independent) 
theory  of  plasticity  formulated  in  a  strain  space  setting  with  the  use  of 
special  constitutive  equations  employed  previously  by  Naghdi  and  Nikkei  (1984). 
Also,  by  employing  the  strain-hardening  characterization  which  arises  in  the 
strain-space  formulation  of  plasticity  (Casey  and  Naghdi,  1981,  1983,  1984a)  and 
which  is  represented  by  a  scalar  function  $>,  we  calculate  (over  the  domain  of 
interest)  the  values  of  <$  for  the  material  used  in  the  experiments  of  Lamba  and 
Sidebottom  (1978a, b)  and  plot  this  as  a  function  of  plastic  strains.  As  in  the 
experiments  of  Lamba  and  Sidebottom  (1978a, b),  our  calculations  are  carried  out 
in  the  context  of  small  deformation. 

*The  abbreviation  OFHC  stands  for  "oxygen-free  high  conductivity." 


In  a  previous  paper  (Naghdi  and  Nikkei,  1984)  comparisons  were  made 
between  the  predictions  of  uniaxial  stress  and  strain  cycling  in  plasticity 
and  corresponding  available  experiments.  The  previous  calculations  (Naghdi 
and  Nikkei,  1984)  were  carried  out  in  a  strain  space  setting  by  means  of 
special  constitutive  equations  obtained  on  the  basis  of  a  theory  first 
formulated  relative  to  stress  space  by  Green  and  Naghdi  (1965,  1966)  and 
subsequently  reformulated  relative  to  strain  space  by  Naghdi  and  Trapp  (1975a), 
along  with  some  additional  developments  pertaining  to  loading  criteria  and 
hardening  characterization  by  Casey  and  Naghdi  (1984a, b). 


Description  of  the  Strain  Cycling  Problem 

As  in  the  experiments  of  Lamba  and  Si  debottom  (1978a,b),we  consider 
the  combined  tension-compression  and  torsion  of  a  thin-walled  circular 
cylinder  whose  axis  is  in  the  x-j- direction.  As  usual,  we  neglect  the  varia¬ 
tion  of  stress  and  strain  in  the  radial  coordinate  direction  so  that  the 
stress  and  strain  components  referred  to  cylindrical  polar  coordinates  are 
then  independent  of  position.  We  denote  the  axial  stress  and  axial  strain, 
respectively,  by  s^  and  e^;  and  similarly  denote  the  shear  stress  and  shear 
strain,  respectively,  by  s^  an<^  e-jg* 

The  analysis  in  Section  3  and  the  procedure  for  calculations  in  Section  4 
require  a  reduction  of  the  general  loading  criteria  since  not  all  components 
of  strain  (and  strain  rate)  are  known  a  priori  in  the  context  of  the 
particular  experiment  considered.  This  reduction  is  discussed  in  Appendix  A 
and  the  results  are  employed  in  all  calculations  of  Section  4.  Our 
calculations  do  not  include  those  appropriate  for  two-dimensional  stress 


cycling,  even  though  it  is  relatively  easier  to  calculate  such  cycles  for 

* 

materials  exhibiting  hardening  behavior. 

In  the  remainder  of  this  section,  we  discuss  in  some  detail  the  main 
features  of  the  calculated  results.  A  summary  of  the  calculations  can  be 
arranged  in  three  groups  corresponding  to  the  nonproportional  strain  cycling 
experiments  of  Lamba  and  Sidebottom  (1978a, b)  as  listed  in  the  opening 
paragraph  of  this  section.  The  first  experiment  pertains  to  saturation  hardening 
after  loading  from  an  undeformed  state,  and  the  second  two  concern  the  behavior 
after  the  material  has  reached  a  state  of  saturation.  In  each  case,  the  strain 
path  in  the  e^-e^  plane  is  the  input  to  the  problem,  and  the  stress  response 
(s^-j  and  s^)  is  calculated  from  the  relevant  constitutive  equations.  A 
graphical  presentation  of  the  calculated  stress  trajectory  in  the  s^-s-^  plane 
alone  does  not  give  all  the  relevant  information,  unless  a  knowledge  of  the 
correspondence  between  all  points  of  the  stress  trajectory  and  all  points  of  the 
input  strain  path  in  the  e^-e^  plane  is  also  known.  It  is  therefore,  necessary 
to  also  plot  the  calculated  results  in  either  the  s^-e^  or  the  2~ei 2  Piflne* 
For  clarity's  sake,  we  have  presented  the  calculated  results  in  both  the  1 

and  si2_ei2  P1anes* 

Saturation  Hardening 

For  two-dimensional  saturation  hardening,  we  prescribe  e^  and  e^ 
parametrically  as  functions  of  time  such  that  in  the  611*612  Plane  the 
strain  path  is  as  depicted  in  the  inset  of  Fig.  1(a).  This 


★ 

This  is  because  for  materials  exhibiting  hardening  behavior  the  loading 
criteria  in  stress  space  may  be  used  (Casey  and  Naghdi,  1984b),  and  for 
two-dimensional  stress  cycling  all  of  the  components  of  stress  are  prescribed 
in  contrast  to  two-dimensional  strain  cycling  where  some  of  the  components 
of  strain  are  not  known  a  priori . 
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corresponds  to  uniaxial  tensile  loading  from  the  undeformed  state  until  the 

axial  strain  reaches  the  value  0.00635  (point  Q  in  Fig.  1(a))  followed  by 

combined  tension-compression  and  torsion  controlled  in  such  a  way  that  e^ 

and  cycle  in  a  counterclockwise  direction  around  the  nearly  elliptical  path  in 

the  inset  of  Fig.  1(a).®  The  strain  components  e^  and  e^  were 

parametrically  specified  to  be  sinusoidal  differing  by  a  90 
degree  phase  shift.  The  calculated  results  for  stresses  are  shown  in 

Figs,  l(a-c)  along  with  the  experimental  curves  for  the  first  few  cycles 
only,  since  the  portion  of  the  trajectory  for  additional  cycles  would  crowd 
the  figure.  Figure  1(a)  represents  the  stress  trajectory  and  Figs.  l(b,c) 
are  alternative  representations  of  the  calculations  in  the  s^-e^  and  $j2~el2 
planes.  Points  corresponding  to  Q,R,$,T  on  the  prescribed  strain  path  in  the 
inset  of  Fig.  1(a)  have  been  labeled  in  Figs,  l(a-c)  only  for  the  portion  of 
the  path  which  represents  a  cycle  at  saturation  and  not  for  prior  cycles. 

These  results  are  in  good  qualitative  agreement  with  the  experimental  results 
of  Lamba  and  Sidebottom  (1978a).  In  agreement  with,  the  experimental  results,  the 
theory  predicts  that  saturation  occurs  after  about  four  cycles.  The  main 
differences  between  the  theoretical  predictions  and  the  experimental  results  are 
that  the  predicted  maximum  shear  stress  at  saturation  (point  R  in  Figs.  1 ( a ,c) ) 
is  slightly  greater  than  the  corresponding  experimental  value  while  the 
predicted  maximum  axial  stress  at  saturation  (point  T  in  Figs.  1 ( a ,b ) )  is 
slightly  less  than  the  corresponding  experimental  value.  The  reasons  for 
these  differences  are  discussed  in  Section  5. 


3 The  actual  experimental  strain  path  in  Fig.  3(a)  of  Lamba  and  Sidebottom  (197ila) 
is  not  perfectly  elliptical.  In  the  present  calculations  the  portion  of  the  path 
in  the  upper  half  of  the  e,,*e12  plane  is  specified  to  be  the  upper  half  or  an 
ellipse  with  vertices  on  the  e,,  axis  at  -0.00623  and  0.00635  and  a  semiaxis  of 
0.00531  in  the  e,-  direction.  The  portion- of  the  path  in  the  lower  half  of  the 
e11-e1?  plane  is  specified  to  be  the  lower  half  of  an  ellipse  with  the  same 
virtiies  on  the  e,,  axis,  but  with  a  semi  axis  of  0.00588  in  the 


:11 


e^  direction. 
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Erasure  of  Memory 

Consider  next  the  erasure  of  memory  property  observed  by  Lamba  and 
Sidebottom  (1978a)  in  the  case  of  OFHC  copper.  The  phenomenon 
occurs  when  the  material  has  reached  a  state  of  saturation 
after  strain  cycling  along  an  elliptical  path  in  the  -j 2  P^ane*  If 
subsequent  linear  strain  paths  lie  within  the  elliptical  strain  path  (which 
was  originally  used  to  reach  the  state  of  saturation),  then  corresponding  to 
a  given  linear  path  in  the  e^-e^  plane  there  is  a  particular  limiting 
stress  cycle  response  in  the  s^-s^  plane.  This  process  is  repeated  as 
the  material  is  cycled  again  along  the  strain  path.  Lamba  and  Sidebottom  (1978a) 

observed  that  a  larger  cycle  (not  necessarily  along  a  linear  path)  in  the 
ell’e12  p^ane  essentially  erases  the  effect  of  any  previously  traversed  smaller 
cycles  and  returns  the  material  to  a  state  in  which  the  limiting  stress  cycle  response 

corresponds  to  that  of  the  larger  cycle.  With  the  use  of  the  linear  strain  paths  BDB  and 
BFB  (as  the  larger  and  smaller  cycles)  indicated  in  the  inset  of  Fig.  2(a), 
a  calculation  leading  to  a  state  of  saturation  was  arrived  at  by  first  cycling 
in  shear  along  B'D^for  four  cycles  followed  by  cycling  along  B'CD'EB'  for  an' 
additional  four  cycles.  Then,  the  effect  of  erasure  of  memory  was  calculated 
by  further  cycling  which  alternates  along  BDB  and  BFB.  The  results  of  calculations 
are  shown  in  Figs.  2(a-c)  along  with  the  experimental  curves  from  Lamba  and 
Sidebottom  (1978a,  Figs.  5(b,c)),  where  for  clarity  only  the  portion  of  the  j 

j 

i 

response  corresponding  to  cycling  which  alternates  along  BDB  and  BFB  are  i 

§§  i 

indicated.  The  erasure  of  memory  phenomenon  is  best  seen  in  Fig.  2(c), where  j 

S  6 

In  Lamba  and  Sidebottom  (1978a),  at  the  point  we’ve  labeled  D,  the  value  of 
the  shear  strain  shown  in  Fig.  (5a)  (the  strain  path  in  the  axial  strain  - 
shear  strain  plane)  is  somewhat  different  from  its  value  shown  in  Fig.  (5b)  1 

(the  experimental  results  shown  in  the  shear  stress-shear  strain  plane).  In  j 

prescribing  the  strain  path  for  our  calculations,  in  order  to  compare  with  I 

the  experimental  results  we  specified  the  value  of  the  shear  strain  at  point  D  ! 

to  be  that  indicated  in  Fig.  (5b)  of  Lamba  and  Sidebottom  (1978a).  1 

i 
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along  the  segment  of  path  BO  the  effect  of  the  path  BFB  has  become  undetectable 
by  the  time  the  trajectory  has  reached  point  D.  Thus,  the  larger  cycle  BOB 
erases  the  material's  memory  of  the  previously  traversed  smaller  cycle  BFB.  In 
addition,  the  overall  response  is  in  qualitative  agreement  with  the  experimental 
results  of  Lamba  and  Sidebottom  (1978a),  while  the  predicted  shear  stress  is 
again  higher  than  tne  corresponding  experimental  results  (Fig.  2(c))  and  the 
predicted  compressive  axial  stress  has  a  larger  absolute  value  than  the 
corresponding  experimental  results  (Fig.  2(b)).  Also,  the  slope  of  the 
calculated  s^2“e-j2  curve  *n  ^9*  2(c)  during  loading  along  F8  is  different 
than  on  the  corresponding  experimental  curve.  Again,  these  differences  are 
discussed  in  some  detail  in  Section  5. 

Complex  Nonproportional  Strain-Path 

Now,  with  reference  to  the  third  experiment  of  Lamba  and  Sidebottom 
(1978b),  consider  the  predicted  response  for  a  "complex"  nonproportional  strain 
path  (using  the  terminology  of  Lamba  and  Sidebottom)  applied  after  the 
material  has  reached  a  state  of  saturation.  The  state  of  saturation  is  attained 
in  the  same  way  as  discussed  earlier  in  the  preceding  paragraph,  namely  by 
cycling  first  in  shear  followed  by  cycling  along  an  elliptical  path  in  the  e^- 
e^  Plane.  After  this,  a  path  was  prescribed  which  returned  e^  and  e^  t0  the 
origin  of  the  e^-e^  P^ane  0  in  Fig.  3(a))  while  at  the  same  time 

the  plastic  strains  were  returned  to  the  value  zero.  Next  let  e^  and  e^  be 


This  calculation  was  performed  by  prescribing  a  path  unloading  from  the  elliptical 
path  at  a  point  such  as  B'  in  the  insert  of  Fig.  2(a)  on  which  e,,=0  and  allowing 
reverse  loading  to  take  place  until  such  a  point  where  e^O.  At  this  point  e?,  has 
also  become  essentially  zero  which  is  again  indicative  of^erasure  of  memory.  Re¬ 
versing  the  direction  along  this  path  until  point  0  in  Fig.  3(a)  is  reached 
results  in  only  elastic  behavior. 
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prescribed  parametrically  as  functions  of  time  such  that  in  the  e^-e^  plane 
the  strain  path  is  represented  by  the  linear  segments  0-1.  1-2,  2-3,  .  .  .,  7-8 
in  the  inset  of  Fig.  3(a).  The  calculated  results  for  stresses  are  shown  in 
Figs.  3(a-c)  which  qualitatively  agree  well  with  the  experimental  results  of 
Lamba  and  Sidebottom  (1978b,  Figs.  2(a-d)).  Moreover,  the  theory  again  predicts 
the  erasure  of  memory  effect  since  as  the  path  0-1,...,  7-8  in  the  inset  of 
Fig.  3(a)  is  traversed  a  second  time,  the  stress  path  becomes  indistinguishable 
from  the  first  time  the  path  was  traversed.  Thus,  the  larger  path  6-7-8  erases 
the  effect  of  the  previously  traversed  smaller  portion  of  the  path  1-6. 

Strain  Hardening  Behavior  -  The  Function  $ 

In  view  of  the  fundamental  role  played  by  a  scalar  function  $  that 

characterizes  strain-hardening  behavior  of  the  material  (Casey  and  Naghdi ,  1984a), 

we  represent  as  a  surface  which  exhibits  the  nature  of  strain-hardening  at  any 

t 

elastic-plastic  state  with  fixed  values  of  the  total  strain  components  e^  and  e12. 

The  surface  in  Fig.  4  is  a  plot  of  4 (defined  in  Eq.  (2.10))  as  a  function  of  the 
plastic  strains  e^  and  e!j2  for  fixed  values  of  the  total  strains  e^  and  e^2 
prescribed  in  the  course  of  calculations;  in  the  plot  of  Fig.  4  the  prescribed 
total  strains  are  specified  to  be  e^=e^2=0.  For  different  values  of  and  e12, 
the  surface  represented  by  merely  translates  parallel  to  the  plane  of  e^-e^. 

A  choice  of  the  total  strains,  say  (e.^,e.|2),  specifies  a  particular  surface 
and  elastic-plastic  states  reached  by  different  strain  paths  ending  at  (e^,e12)  will 
in  general  correspond  to  different  points  on  the  surface  4.  The  outer  boundary  of  this 
surface  (which  has  the  smallest  valueof$)  corresponds  to  the  state  of  saturation,  whil 
the  inner  boundary  corresponds  to  the  largest  value  that  i>  can  take  in  the  domain 

$A  topographical  representation  of  this  kind  was  used  recently  in  a  different 
context  by  Casey  and  Lin  (1983). 
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of  interest.  A  further  discussion  of  how  the  surface  was  calculated  is  found  at 
the  end  of  Section  4.  As  is  evident  from  the  plot  in  Fig.  4,  <t>  is  always 
positive;  and  hence,  in  view  of  the  conditions  for  strain-hardening  characterization 
(see  for  example  the  conditions  (8)  in  Naghdi  and  Nikkei ,  1984),  the  material 
always  exhibits  hardening  behavior. +  It  also  shows  how  4  decreases  with  addition?! 
plastic  deformation,  taking  its  largest  value  at  initial  yield  and  its  smallest 
value  at  saturation.  The  value  of  $  is  constant  at  both  initial  yield  and 
saturation. 


+It  should  be  emphasized  that  even  though  one  component  of  the  stress  may  be 
decreasing  during  loading  while  the  corresponding  strain  component  is 
increasing  (e.g.,  as  in  Fig.  2(b))  the  material  is  not  exhibiting  softening 
behavior  in  two-dimensional  cycling. 
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With  reference  to  a  strain-space  formulation  of  plasticity  and  for 
the  special  constitutive  equations  used  previously  (Naghdi  and  Nikkei  1984), 


we  include  here  a  brief  summary  of  the  relevant  equations  of  the  purely 
mechanical  theory  contained  in  the  papers  of  Green  and  Naghdi  (1965,  1966), 
Naghdi  and  Trapp  (1975a)  and  Casey  and  Naghdi  (1981,  1983)'. 

In  the  context  of  infinitesimal  deformation,  we  recall  that  the 


main  ingredients  of  the  rate-type  theory  of  plasticity,  in  addition  to  the 
total  strain  e^,  are  plastic  strain  ej^  and  a  measure  of  work-hardening  k. 
Also,  no  distinction  needs  to  be  made  between  various  measures  of  stress 
which  we  denote  by  s^.  As  usual,  it  is  convenient  to  express  the  various 
constitutive  response  functions  in  terms  of  the  components  of  deviatoric 


stress  and  deviatoric  strain 


Y 


KL 


,  namely 


KL  =  SKL  ' 

s  6kl 

’  tkk 

KL  =  eKL  " 

iSKL 

*  ykk 

YKL  eKL  "  e  5KL  ’  YKK  =  0  ’ 

with  a  similar  definition  for  the  deviatoric  plastic  strain  yj^,  where 
?,  e  and  e^  denote  the  mean  normal  stress,  mean  normal  strain,  and  mean 
normal  plastic  strain,  respectively.  For  materials  which  are  isotropic 


in  reference  configuration  and  in  the  presence  of  plastic  incompressibility 
(I*5  =  0),  we  specify  the  stress  response  by  generalized  Hooke's  law  and  the 
loading  functions  g  and  f  in  strain  space  and  stress  space,  respecitvely ,  by ® 


TA  more  expanded  summary  is  contained  in  Sections  2  and  3  of  Naghdi  and  Nikkei 
(1984). 

§ 

The  loading  functions  in(2.2)with  a  being  regarded  as  constant  (rather  than  a 
function  of  <)  were  used  in  the  paper  of  Caulk  and  Naghdi  (1978)  and,  with  o 
as  a  function  of  k,  are  the  same  as  those  used  by  Naghdi  and  Nikkei  (1984). 


g  -  4u2[ykl-(i  +4^-)y5i  Ky,,  -0  +*#)-£,  ] 


4m  ;,KLjl,KL 


4 u  ' -KL- 


^TKL~2  ^<^yKL^tKL’2  5(<)yKL}  "  K  “  f 


(2.2) 


with 


(a.  -  a. )<  +  a.ic  -  a  <_ 
a  •  a(<)  «  -2 - S - s-° - 0  s 


Ko  ‘  *s 


(2.3) 


In(2.2)  and  (2.3)  the  coefficients  aQ,  a,  .  <Q,  are  constants  and  u  is  the 

elastic  shear  modulus.  The  constants  ac  and  as  are  so  chosen  that  ot{<)  takes 

the  value  a  when  <  =  <„  and  takes  the  value  cr  when  <  -  k„. 

oo  s  s 

We  adopt  the  loading  criteria  or  strain  space  as  primary.  Then, 

after  invoking  the  work  assumption  o-r  Naghdi  and  Trapp  (1975b),  the 

constitutive  equations  for  the  rate  ■'if  plastic  strain  e£L  and  the  rate  of 

work- hardening  parameter  <  may  be  expressed  as  (see  the  development  among 

Eqs.  (36) -(42)  in  Casey  and  Naghdi  1981): 

0,  when  g  <  0,  or  when  g  =  0  and  g  <  0, 


'KL 


rP  - 
rKL 


-ayPL)  ,  when  g  =  0  and  g  >  0 


(2.4a,b) 


and 


<  s  (B(<)tkl  +  h(<h^L)ePL  . 


In  (2.4)  and  (2.5),  the  quantities  B  '  are  defined  by 


(2.5) 


8  =  B(<)  *  (^-)6 

Ko  s 


n  *  n(<)  =  ( 


K 

0  s 


)n 


(2.6) 


B,n  are  constants. 


s_la_ 


9  *  ^®KL  *KL  =  4U^TX1-"  2°r‘  L^KL 


(2.7) 


and  on  the  yield  surface  g  *  C, 
A.  =  8ptc  >  0  , 

0  S  i 


r  ■  2SKt2[1t(^)l%-ti'KL)4](2THN^-4)(TW-|4) .  (2.9) 


(2.8) 


0  s 
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The  yield  function  f  in  (2.2)  is  of  the  von  Mises  type:  it  is  quadratic 
in  the  deviatoric  stresses  and  allows  for  translation  and  change  in  size  of 


the  yield  surface.  The  amount  of  translation  of  the  yield  surface  is  linear  in 
the  deviatoric  plastic  strain  and  is  also  linear  in  <  through  the  coefficient 
3  defined  by  (2.3).  Also,  the  coefficient  functions  §(<)  and  fj(ic)  which  occur  in 
(2.5)  are  linear  in  <,  they  assume  the  respective  values  6  and  n  at  initial  yield 
when  <  =  <Q,  and  both  vanish  at  saturation,  i.e.,  when  c  *  <s>  <  =  0. 

We  also  recall  here  that  the  strain-hardening  behavior  may  be  characterized 

by  means  of  a  rate- independent  scalar  quantity  denoted  here  by  $  (Casey  and 

*  a  ^ 

Naghdi ,  1981,  1983).  Ouring  loading  the  quantity  $  has  the  same  value  as  f/g, 

where  ?  5  (af/Ss^s^  and  $  is  defined  by  (2.7)^.  With  the  particular 
constitutive  equations  (i.e.  generalized  Hooke's  law)  used  here  and  the  work 
assumption  of  Naghdi  and  Trapp  (1975b),  $  can  be  expressed  as 

(2.10) 

and  in  fact  since  r  +A  >0,  r  alone  may  be  used  to  characterize  the 
strain-hardening  behavior  (Casey  and  Naghdi  1984b,  Eqs.  (4.50)  and  (4.51)). 


Tt  was  demonstrated  recently  by  Casey  and  Naghdi  (1984a,  Eq.  (3.8))  that 
the  quantity  $  is  equal  to  the  determinant  of  a  certain  rate-independent 
fourth-order  tensor  which  plays  c  fundamental  role  in  the  theory  of  plasticity 
and  which  arises  naturally  in  relating  the  time  rate  of  stress  to  the  time  rate 
of  strain. 


Equations  for  a  Two-Dimensional  State  of  Stress 


Consider  now  the  homogeneous  deformation  of  an  elastic-plastic 
material  sustained  by  a  biaxial  state  of  stress  (one  normal  component  and 
one  shearing  component)  in  which  the  corresponding  two  components  of  strain 
are  prescribed  functions  of  time.  Lamba  and  Sidebottom  (1978a, b)  modelled 
this  experimentally  by  performing  a  strain-controlled  combined  tension- 
torsion  test  of  a  thin-walled  circular  cylinder  in  which  the  axial  strain 
and  the  shear  strain  were  prescribed  functions  of  time.  With  the  notation 
s11(t)  for  the  axial  stress  and  s-^t)  for  the  shear  stress,  the  two- 

dimensional  state  of  stress  can  be  represented  in  matrix  notation  as 

NTKlJI  a  ^ll  +s12  H  aKL  N  ’  s  =  3*11 

where  the  constant  matrices  a„.  and  b..  defined  by 


.  I  0  1  0  |  z  v  u  I 

aK,il  ■  |  1  0  0  ,  Ijb,..  ||  *  0-10 

L  1  0  0  0  1  KL  0  0-1  1 


2  0  0 


are  introduced  for  convenience. 

The  intended  calculations  require  prescribing  e^  and  e^  parametrically 
as  particular  functions  of  time  (which  specify  a  path  in  the  e^-e^  Plane), 
while  the  other  two  components  of  normal  strain  and  833)  remain  so  far 
unspecified  and  hence  unknowns.  Since  ej^  *  633  and  8^3  =  ®23  3  0 

by  (2.40  and  (3.1),  it  follows  that  the  plastic  strain  tensor  can  be  expressed 


in  the  form 


11 


II  £  11 


=  Iep 
V 11 


where  ep-j  and  e^  depend  only  on  time.  From  the  stress  response  (generalized 
Hooke's  law)  we  have 


s. 


111 


where  £  is  Young's  modulus.  From  the  inverted  form  of  the  stress  response. 


as  well  as  (3.$  and  (3.4)^,  and  after  adopting  the  notation  e„  =  e 


22  *  e33 


we  obtain 


en  ei2 0 

e12  er  0 
0  0  e_ 


(3.5) 


where 


er  *  "I  S11  '  2e11  *'vell  “2^~2v)eii  (-3.6) 

and  where  v  s  (E-2u)/2 y  is  Poisson's  ratio.  Equation  (3.6)  relates  the 
unspecified  components  of  strain  to  the  axial  strain  (which  is  known)  and  the 
axial  component  of  plastic  strain.  But  ef  is  still  unknown,  inasmuch  as 
remains  unknown  until  all  of  the  constitutive  equations  pertaining  to  a 
particular  two-dimensional  strain  path  have  been  integrated. 

In  order  to  simplify  some  of  the  expressions  that  follow,  we  introduce 
the  abbreviations 


M  =  s.j  i  -  ^  a(<)  e^-j  »  N  =  si2  ’J“(<)  e12 


(3.7) 


Then,  with  the  use  of  (3-1)  to  (3.6),  the  loading  functions  in(2.2)now  assume 
the  simplified  forms 


9  .|E2[e,r(lt||)e?,]zt8uz[e12-(lt|;)eP2]z-.  . 


3» 

a 


f  -  fn2* 


2N2  -  < 


(3.8) 


For  the  two-dimensional  state  of  stress  defined  by  (3.1),  A  is  still  given 
by  (2.8);  but,  on  the  yield  surface,  r  in  (2.9)now  assumes  the  special  form 

r  «2a<+3[l  +  (^5-XMeP1  +  2Ne52)][M(Ss11+mi)+  3N(8s12  +  ne^2)].  (3.9) 
Ko  s 

In  view  of  (3.1)  ,  (3.3),  and  (3.5),  from  (2.7)  the  quantity  g  occurring  in  the 
loading  criteria  in  strain  space  is  given  by 
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(3.10a 


=  +  8y  +  y(3y-E)Me^  , 

and  the  corresponding  quantity  f  in  stress  space  is 


f  =  |msh  +  4Ns]2  #  (3.11) 

where  in  obtaining (3. 10b) the  time  derivative  of  (3.6)  has  been  used.  It 
should  be  noted  that  g  depends  on  the  unspecified  strain  rate  ep  (or  the 
unknown  axial  plastic  strain  rate  4^ ) ;  and  hence  the  determination  of 
whether  the  state  of  strain  at  an  elastic-plastic  state  is  undergoing  loading, 
neutral  loading,  or  unloading  cannot  in  general  be  ascertained  directly  from 
(3.10). 

Since  during  loading  the  constitutive  equation (2. 4b)  for  the  plastic 
strain  rate  c'jnt-iins  g  on  its  right-hand  side,  in  view  of  (3. 10b)  it  is  at 
once  obvious  that  e^  occurs  on  both  sides  of  the  equation  in  the  case  of 
axial  plastic  strain  rate.  After  solving  this  equation,  a  new  form  of  the 
constitutive  equation  results  which  depends  only  on  the  specified  strain 
rates  e^  and  e^.  This  equation  then  enables  the  constitutive  equations 
for  the  other  components  of  plastic  strain  rate  to  be  expressed  in  terms  of 
only  e^  and  e-^  also.  This  procedure  along  with(3.1)  and  (3.3)  leads  to  the 
following  expressions  for  the  non-zero  plastic  strain  rates  during  loading 
(g-0,  g  >  0) 


•p  _  4  M 
"11  ”  3  Q  9 


ip  =ap  = 
e22  e33 


15P 

2  el  1 


where  we  have  introduced  the  quantities  Q  and  g  defined  by 


(3.12) 


Q  s  r  +  A  -  t£(3u-  e)m2 


(3.13) 


It  is  shown  in  Appendix  A  under  the  assumption  that  Q  >  0  always,  which 
represents  a  range  of  strain-hardening  behavior  sufficiently  general  for 
our  present  purposes,  that  a  knowledge  of  only  the  prescribed  strain 
rates  e^  and  e^  is  sufficient  to  determine  whether  the  material  at  an 
elastic-plastic  state  is  undergoing  loading,  neutral  loading,  or  unloading. 
In  fact,  with  Q  >  0  the  quantity  g  may  be  used  in  special  loading  criteria 
appropriate  for  the  particular  problem  under  discussion.  Thus,  with  Q  >  0, 
the  constitutive  equations  for  the  plastic  strain  rate  and  the  rate  of  work 
hardening  may  alternatively  be  expressed  as+ 

0  ,  if  g  <  0,  (a) 

0  ,  if  g  =  0  and  g  <  0,  (b) 

•P  = 

eKL  0  ,  if  g  *  0  and  §  =  0,  (c) 

Eq. (3.72) ,  if  g  =  0  and  3  >  0,  (d) 

and 

K  *  (Bs^  +|  ne^  )e^  +  2(§s.j2  +  n  e^2)  e^2  , 

where  the  loading  criteria  in  (3.15a, b,c,d)  correspond,  respectively,  to  an 
elastic  state,  unloading  from  an  elastic-plastic  state,  neutral  loading  at  an 
elastic-plastic  state,  and  loading  at  an  elastic-elastic  state. 

We  note  that  during  unloading  and  neutral  loading  it  follows  from  (2.4a), 
(3.4),  (3.10b)  and  (3.11)  that  g  =  §  =  ?,  while  during  loading  from  (3.4) 
and  (3.11)  to  (3.14)  we  have 

f  =  y  3  *  7~irK  9  *  ^9=0  »  9  >  °^  • 

TWe  emphasize  that  the  loading  criteria  in  (3.1$  derived  from  the  general 
form  of  the  strain  space  loading  criteria,  and  not  postulated  in  an  ad 
hoc  fashion. 


Results  from  the  general  theory  similar  to  these  and  those  obtained  in 
Appendix  A  are  summarized  in  Table  1  of  Casey  and  Naghdi  (1984c). 

8efore  closing  this  section,  we  need  to  corrment  further  on  the 

restriction  Q  >  0  which  bears  on  the  loading  criteria  in  (3.19.  With 

reference  to  (3.13),  since  r  may  be  used  to  characterize  the  strain-hardening 

behavior, with £.£}  and(3.8)  the  condition  Q  >  0  may  be  seen  to  include  hardening 

behavior  (<5>0)  and  perfectly  plastic  behavior  ($*0),  but  at  first  sight  may 
appear  to  exclude  a  small  range  of  softening  behavior.5  However,  it  can 

be  shown  (Appendix  A)  that  at  an  elastic-plastic  state  with  Q  <  0,  the 

applied  strain  rates  e^  and  e^  can  only  be  such  that  g  <  0.  But  the 

expression  (3.14)2  which  occurs  in  the  loading  criteria  (3. 15) must  be  capable 

of  admitting  all  possible  choices  of  and  §12  (  and  hence  capable  of 

taking  both  positive  and  negative  values)  for  all  physically  realistic 

tests.  Since  a  state  corresponding  to  Q  <  0  is  unnecessarily  restrictive, 

it  will  be  excluded  from  consideration  in  the  present  development. 


This  is  because  the  condition  Q  >  0  places  a  greater  restriction  on  r  than 
does  the  more  general  result  r  +  A  >  0  (which  is  a  consequence  of  the  work 
assumption  of  Naghdi  and  Trapp,  1975b). 
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4.  Determination  of  Material  Constants-  Details  of  Calculations. 

Previously  a  procedure  was  suggested  for  determining  the  material 
coefficients  in  the  constitutive  equations  from  the  experimental  data 
in  a  uniaxial  cyclic  loading  test  (Naghdi  and  Nikkei,  1984).  This 
procedure  is  by  no  means  unique;  and,  in  fact,  for  data  obtained  from  one¬ 
dimensional  tests,  other  than  uniaxial  cyclic  tests,  it  may  be  desirable 
to  use  a  different  procedure.  Such  alternative  procedures  are  likely  to  be 
more  advantageous  in  two  or  three  dimensional  tests.  With  reference  to 
their  experimental  results  for  OFHC  copper,  Lamba  and  Sidebottom  (1978a) 
note  that  the  peak  axial  stress  attained  after  the  material  has  been  cycled 
along  a  strain  path  such  as  that  in  Fig.  1(a)  is  significantly  higher  than 
the  peak  axial  stress  attained  after  simple  uniaxial  cycling.  They  go  on  to 
state  that  this  difference  indicates  that  ".  .  .  material  properties 
obtained  from  tensile  tests  or  even  uniaxial  cyclic  tests  will  not  give 

accurate  predictions  of  cyclic  deformation  under  nonproportional  or  out-of- 
★ 

phase  conditions."  A  conclusion  of  this  kind  cannot  be  made  independent  of 
a  particular  theory  used  and  does  not  follow  from  the  knowledge  of  only  a 
part  of  the  experimental  measurements  (e.g.,  a  part  of  the  measurements  that 
leads  to  the  calculated  s^).  Rather,  it  requires  a  detailed  examination  of 
the  entire  experimental  data  in  conjunction  with’ the  relevant  constitutive 
equations. 

As  can  be  seen  from(3.7)-|  2  and(3.8)2  with  <  =  <$ ,o(k)  s  ag  and  with  f=0, 
the  expression  for  in  a  two-dimensional  cycling  test  depends  on  s^,  s^* 
e^,  e1^  as  well  as  a*a$.  However,  the  expression  for  in  the  case  of  a 
a  "peak"  axial  stress  (i.e.,  the  value  of  s^  when  s12  *  0)  will  depend  only 
on  the  axial  stress  s^  in  addition  to  other  parameters  such  as  a$  and  the 

★  -  i  - 

The  term  out-of-phase  refers  to  the  specification  of  e,,  and  e,2  as  sinusoidal 
functions  of  time  with  the  90  degree  phase  difference  between  them,  resulting 
in  an  elliptic  path  in  the  e^-e^  Plane. 
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plastic  strains.  We  further  observe  that  in  the  biaxial  tests  under  discussion 
the  value  of  the  parameter  <s,  which  represents  the  size  of  the  loading  surface 
in  two-dimensions,  is  particularly  sensitive  to  the  values  of  the  stresses  but 
much  less  so  to  the  values  of  as  and  the  plastic  strains.  In  view  of  the 
experimental  result  of  Lamba  and  Sidebottom  (1978a)  mentioned  in  the  previous 
paragraph,  if  is  determined  from  uniaxial  cycling  data  alone,  then  the 
theoretical  prediction  for  the  two-dimensional  cycling  experiments  may  not  be 
as  accurate  as  when  <s  is  determined  from  the  entire  experimental  data  of  a 
two-dimensional  strain  cycling  test. 

We  keep  the  foregoing  discussion  in  mind  when  describing  the  alternative 
procedure  for  determining  the  material  coefficients  in  the  constitutive 
equations.  A  value  for  <s  will  be  determined  from  two-dimensional  strain 
cycling  experiments.  For  our  present  purpose,  it  will  suffice  to  determine 
the  other  constants  from  uniaxial  cycling  data  as  before.7  Thus,  a  value  for 
<Q  follows  from  f=0  and  the  yield  stress  in  uniaxial  tension,  and  a  value 
for  a$  can  be  determined  from  the  slope  of  the  uniaxial  stress  strain  curve 
at  saturation.  By  selecting  a  point  at  saturation  from  the  experimental  data 
for  a  two-dimensional  cycling  experiment,  similar  to  points  Q,  R,  S  or  T 
shown  in  the  calculated  results  in  Figs.  1  (a-c )  'and  with  the  values  of  stress 
and  plastic  strain  at  such  a  point  of  value  for  <$  may  be  determined  from  f*0. 

We  took  the  value  <s  to  be  the  average  of  the  values  obtained  from  points  R  and  T, 

in  order  to  have  the  best  overall  prediction  for  all  loading  paths  in  the  e^-e^  plane. 
The  remaining  constants  may  then  be  determined  from  the  uniaxial  cycling  data 

by  the  procedure  described  previously  (Naghdi  and  Nikkei,  1984).  Using  the 
above  procedure  with  the  experimental  data  of  Lamba  and  Sidebottom  (1978a, 

Figs.  2  and  3(a-d))  for  OFHC  copper  the  material  constants  were  determined  to 

fIn  contrast  to  for  the  relatively  simple  constitutive  equations  used  here 
the  determination  of  the  other  constants  from  the  uniaxial  tests  seems  to  be 
adequate. 
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be 

J-  =  y-  =  8.14*  10"3  ,  |=2.30*10'2  ,  -9.38  *  10"5  , 

K  tC 

-|=  4.03 *  TO"8  ,  -|=  1.93  *10'6  ,  |=0.374  ,  (4.1) 

E  =  115  GPa  . 

In  the  course  of  identifying  values  for  the  coefficients  aQ  and  ct$,  it  was 
found  that  their  values  differed  only  by  less  than  0.3%.  This  suggests  that 
if  this  difference  can  be  neglected,  we  could  set  oQsas  approximately  and 
then  the  coefficient  n  in (2.5)  and  6) will  no  longer  require  an  independent 
identification.  This  would  make  the  task  of  determination  of  the  coefficients 
much  easier;  and,  in  fact,  with  aQ=as  the  coefficient  n  can  be  determined  in 
terms  of  aQ  and  6.  To  see  this,  we  recall  that  in  the  special  case  in  which 
aQ  =  as  =  a  (say)  the  constitutive  equations  used  here  reduce  to  a  special  case 
of  those  employed  previously  by  Caulk  and  Naghdi  (1978,  Eqs.  (40),  (56)  and  (70^). 
For  their  equations,  they  obtained  the  restriction  n=-|a8  (see  Eq.  (7$,  in  Caulk  and 
Naghdi,  1978).  In  view  of  the  fact  that  the  values  determined  for  aQ  and  ot$  are  essen¬ 
tially  the  same,  for  the  purpose  of  comparison  with  the  experimental  data  of  lamba  and 

Si  debottom  ( 1978a ,  b)  for  OFHC  copper,  it  will  suffice  to  take  a  =  ct  and 

o  s 

n  =  -  ^o6’ 

The  calculations  summarized  in  Section  1  were  carried  out  by  first 
parametrically  specifying  the  strains  e^  and  e^  as  functions  of  time, 
corresponding  to  a  particular  path  in  the  e^-e^  Plane.  The  constitutive 
equations  (3.15)  and  (3.16)  were  then  integrated  numerically  with  the  values  (4.1) 
for  the  coefficients.  Finally,  the  stresses  were  calculated  using  (3.4). 

To  calculate  <f  in  (2.10),  a  knowledge  of  A  in  (2.0)  and  F  in  (3.9)  is 
needed.  The  particular  representation  of  r  in  (3.9)  depends  on  the  stresses. 
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as  well  as  the  quantities  M  and  N  which  also  involve  the  stresses  s^  and  s^. 
Then,  in  view  of  (3.4),  r  may  be  expressed  in  terms  of  the  total  and  plastic 
strains.  It  follows  that  4  depends  on  total  and  plastic  strains,  as  well  as  on  < 
But,  from  the  fact  that  on  the  yield  surface  g  =  0  with  the  function  g  given  by 
(3.8)-],  it  is  possible  to  solve  for  <  in  terms  of  the  other  quantities  so  that 
$  may  be  represented  as  a  function  of  only  e^,  e.^,  e^,  and  e^2*  It  is  of 
interest  to  plot  the  variation  of  4  with  plastic  strains  for  fixed  values  of 
e^-j  and  e^.  After  making  the  substitutions  indicated  above,  the  dependence 
of  4  on  the  total  and  plastic  strains  may  be  expressed  in  the  form 

4  1+J^)'1*,,  .  e\2  -  (4.2 

For  definiteness,  we  specify  e^  a  e12  =  0  and  then  calculate  the  value  of 
4  for  each  e^,  e^2  pair.  The  resulting  plot  using  the  values  (4.1)  is 
shown  in  Fig.  4.  It  is  clear  from  the  arguments  of  the  function  4  in 
(4.2)  that  this  single  calculation  provides  all  of  the  relevant  information 
on  the  variation  of  4;  and,  for  any  other  specified  values  of  e^  and  e12> 
the  surface  plotted  in  Fig.  4  will  not  change  in  shape  but  will  simply 
translate  parallel  to  the  e^  -  e^2  plane  by  the.  constant  amounts 

<'+TT>'\l  •  (,*^>'le12  <4'3 

in  the  e^  and  e^2  directions  respectively.  It  must  be  kept  in  mind  that 
the  domain  in  the  e^-e^2  P^ane  must  be  such  that  the  value  of  <  corresponding 
to  the  values  of  e^,  e^2,  e^  and  e^2  is  between  <Q  and  <$.  In  Fig.  4  the 
outer  boundary,  with  the  smallest  value  of  4  corresponds  to  k  *  <s  and  the 
inner  boundary,  with  the  largest  value  of  4  corresponds  to  <  *  kq.  It  should  be 
mentioned  that  some  elastic-plastic  states  corresponding  to  points  on  the  surface 


may  not  be  reached  by  any  path.  For  example,  on  the  inner  boundary  in  Fig.  4, 
where  <  *  <Q,  the  plastic  strain  is  nonzero  indicating  some  plastic  deformation 
must  occur  to  reach  an  elastic-plastic  state  corresponding  to  this  edge  of  the 
surface.  However,  after  plastic  deformation  has  taken  place  <  cannot  have  the 
value  <q  (the  vale  of  <  at  initial  yield)  while  the  material  exhibits  hardening 
behavior  except  in  the  special  case  of  purely  kinematic  hardening.  Thus,  for  the 
plot  displayed  in  Fig.  4,  no  actual  elastic-plastic  state  corresponding  to  the 
inner  boundary  of  the  surface  can  be  reached. 


5.  Concluding  Remarks 

In  summary,  a  relatively  simple  set  of  constitutive  eouations  is  used  to 
predict  various  phenomena  occurring  in  two-dimensional  strain  cycling  in  the 
range  of  small  deformation.  The  calculated  predictions  are  compared  with 
corresponding  experimental  results  of  Lamfaa  and  Sidebottom  (1978a, b)  with 
good  qualitative  agreements.  It  is  noteworthy  that  even  during  the  post-saturation 
behaviour  of  the  material,  the  constitutive  equations  used  have  adequate  predictive 
capabilities  as  demonstrated  by  the  results  in  Figs,  (2(a-c))  and  Figs.  (3(a-c)). 
Also,  it  may  be  emphasized  that  the  theoretical  calculations  successfully  predict 
the  erasure  of  memory  phenomenon  which  has  significant  practical  utility  in  two- 
dimensional,  post-saturation  strain  cycling  experiments,  dearly,  by  exploiting 
this  phenomenon,  the  response  of  a  material  to  strain  cycling  along  strain  paths 
in  different  directions  can  be  determined  by  performing  experiments  on  one  specimen 
if  a  large  cycle  which  erases  the  material's  memory  of  the  preceding  smaller  cycle 
is  traversed  after  each  smaller  cycle. 


The  main  differences  between  the  present  calculations  and  the  experimental 
results  can  be  attributed  primarily  to  the  manner  of  identification  of  two  of  the 
material  constants  at  saturation,  namely  <s  and  a  ,  These  values  can  be  chosen 
to  be  the  same  as  those  identified  either  from  experiments  in  simple  tension 
(e^  3  0)  or  experiments  in  simple  shear  (e^  =  0),  .  The  values  resulting  from 
experiments  in  simple  tension  can  lead  to  discrepancies  in  matching  experiments 
involving  mainly  simple  shear  (such  as  those  in  Figs.  (2 (a-c))  and  (3(a-c)),  and 
likewise  values  chosen  from  experiments  in  simple  shear  will  affect  the  agreements 
between  the  theoretical  and  experimental  results  in  simple  tension.  The  choice  for 
the  value  of  the  constant  used  in  the  calculations  (for  details  see  Sec.  4) 
was  motivated  by  a  desire  to  obtain  the  best  overall  agreement  with  experimental 
results  for  all  directions  of  the  strain  paths  in  the  e^-e^  plane.  Similarly, 
as  noted  in  the  last  section,  a$  was  chosen  to  be  equal  to  aQ  in  order  to  simplify 


the  identification  of  the  material  coefficients  for  OFHC  copper.  The  slope  of  the 


calculated  response  curve  in  the  s-\2~e]2  P^ane  at  saturation  can  be  shown  to 

depend  to  some  degree  on  a$,  so  that  if  a  more  accurate  prediction  of  the  slope  of 

the  curve  during  loading  along  FB  in  Fig.  2(c)  is  desired  an  alternative  procedure 

★ 

for  determination  of  a$  could  be  used. 

The  plot  displayed  in  Fig.  4  contains  detailed  information  pertaining  to 
strain-hardening  behaviour  for  the  particular  material  used  in  the  experiments 
of  Lamba  and  Sidebottom  (1978a, bj.  Accessibility  to  such  information  or  the  data 
representing  $  for  all  values  of  total  and  plastic  strains  in  the  domain  of 
interest  is  clearly  of  value  in  analyses  and  computations.  It  should  be  possible 
in  principle  to  determine  the  value  of  $  directly  from  experiments.  Casey  and 
Naghdi  (1984b,  Eq.  (4.29))  have  previously  shown  that  the  function  <t>  can  be 
interpreted  in  terms  of  the  ratio  of  the  outward  normal  velocities  of  the  yield 
surfaces  in  stress  space  and  in  strain  space.  It  is  suggested  that  future 
experimenters  provide  direct  measurements  concerning  the  yield  surfaces  in  both 
stress  space  and  strain  space  and  also  consider  the  possibility  of  obtaining  the 

values  of  $  directly  in  the  course  of  their  experiments. 

We  close  this  section  with  some  remarks  concerning  additional  experimental 
data  on  two-dimensional  cycling  that  have  become  available  very  recently. 


In  addition  to  the  two  elastic  constants,  there  ara  six  material  constants 

U0>as>3>n»<0»<s)  which  must  be  determined  from  experimental  results.  However, 

for  the  special  case  in  which  a  =  a  as  utilized  in  Section  4,  an  additional 

o  s 

restriction  obtained  by  Caulk  and  Naghdi  (1978)  can  also  be  invoked  and  this 
reduces  the  number  of  independent  constants  to  be  determined  to  the  four 
(a0’8>K0*<s) •  With  a0  f  the  calculations  will  involve  more  complex 
expressions  and  plots  (for  different  values  of  e^.e^lof  the  function  $  will 
exhibit  changes  in  shape  and  will  not  simply  represent  translations  of  the 
surface  parallel  to  the  e^-e^  plane  as  in  the  case  when  aQ  -  a$. 


McDowell  (1985)  has  reported  experimental  results  for  two-dimensional  strain  cycling 
of  a  type  similar  to  those  of  Lamba  and  Sidebottom  (1978a, fc)’,  and  in  principle, 
similar  comparisons  can  be  made  with  his  data.  However,  the  data  provided  in  McDowell's 

paper  are  insufficient  for  the  identification  of  all  of  the  material  constants 

s 

appearing  in  the  constitutive  equations  used  here.  Also,  it  should  be  noted  that  our 
approach  for  theoretical  predictions  differs  from  the  "two  surface  stress  space 
model"  used  by  McDowell  (1985)  for  comparison  with  his  experiments, 

Ohashi  et  aJL  (1985)  have  reported  some  experimental  data  for  two-dimensional 
stress  cycling.  They  do  not  discuss  any  theoretical  predictions  for  comparison 
with  their  experiments.  While  in  principle,  there  should  be  no  difficulty  in  making 
such  comparisons,  the  data  provided  is  again  insufficient  for  the  identification  of 
all  of  the  material  constants  appearing  in  the  constitutive  equations  used  here.* 
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^McDowell's  paper  (1985)  does  not  include  the  data  for  uniaxial  strain  cycling  needed 
for  the  identification  of  material  constants  in  the  context  of  the  procedure  used  here. 

+0hashi  et  aK  (1985)  do  not  perform  any  uniaxial  strain  cycling  tests.  The  data 
from  such  tests  are  necessary  for  identification  of  most  of  the  material  constants 
using  the  procedure  s unman  zed  in  Sec,  4. 
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Appendix  A 

We  include  here  the  details  of  the  arguments  which  demonstrate  that 
the  special  loading  criteria  in(3.l5)with  g  (rather  than  g)  are  consistent 
with  the  loading  criteria  of  the  strain-space  formulation,  i.e.,  at  an 
elastic-plastic  state  with  Q  >  0  we  establish  the  correspondence 

g  <  0  •  g  <  0  ,  (unloading), 

g  =  0  •  g  =  0  ,  (neutral  loading), 

g  >  0  •  g  >  0  ,  (loading)  . 

We  first  prove  that  at  an  elastic-plastic  state  with  Q  f  0,  g  *  0 
if  and  only  if  g  =  0.  Sufficiency  follows  from  (3. 12)  .where  if  §  *  0  so  is 
ef,  =  0  and  then  from  (3.14)  g  =  0  also.  To  establish  necessity  suppose 
that  g  *  0.  It  then  follows  from  (2.4)  and  (3. 14}  that  g  =  0. 

In  view  of  the  result  of  the  previous  paragraph  note  that  at  an 
elastic-plastic  state  with  Q  t  0,  g  t  0  if  and  only  if  g  t  0.  We  now 
prove  that  at  an  elastic-plastic  state  with  Q  >  0,  §  >  0  if  and  only  if 

g  >  0.  To  establish  sufficienty  suppose  that  g  >  0.  If  g  <  0,  then  eJJL=  0 

by  (2.4)  and  (3.14)  implies  §  =  g,  which  is  a  contradiction  and  hence 
we  must  have  g  >  0.  To  establish  necessity,  we-suppose  that  g  >  0,  and 
this  implies  that  g  =  +A)/Q  from  (3.17").  In  view  of  the  fact  that 

r  +  A  >  0  (Casey  and  Naghdi,  1984b,  Eq.  (4.50))  and  the  fact  that  we  are 

considering  only  the  case  in  which  Q  >  0  this  implies  g  >  0. 

As  a  consequence  of  the  results  of  the  preceding  two 
paragraphs,  it  follows  that  at  an  elastic-plastic  state  with  Q  >  0,  §  <  0 
if  and  only  if  g  <  0. 

In  the  remainder  of  this  appendix  we  demonstrate  that  if  Q  in  (3.13)  is 
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nonpositive  at  an  elastic-plastic  state,  i.e.,  if  Q  <  0,  then  the  applied 


strain  rates  and  e^  can  only  be  such  that  g  <  0. 

We  first  prove  that  with  Q  <  0  we  can  only  have  g  <  0.  Suppose  g  >  0. 
If  g  <  0,  then£-4)  and  (3.14)  imply  g  =  g  which  is  a  contradiction. 
Alternatively  if  g  >  0,  then  g  a  g(r  +  A)/Q  from  (3.17).  In  view  of  the 
fact  that  r  +  A  >  0  (Casey  and  Naghdi,  1984b,  Eq.  (4.50))  and  that  we  are 
considering  the  case  in  which  Q  <  0,  it  follows  that  £  <  0  which  is  also  a 
contradiction.  Hence,  we  can  only  have  g<0if  Q  <  0. 

Next,  we  prove  that  ifQ=0,  we  can  only  have  g  <  0.  Again,  suppose  that 
g  >  0.  If  g  <  0,  then  (2.4)  and  (3.14)  imply  g  =  g  which  is  a  contradiction. 
Alternatively,  if  g  >  0,  with  the  use  of  the  expression  for  g  resulting  from 
(3 . 1 4)i and  the  identity  resulting  from(3.1$  after  setting  its  left-hand  side 
equal  to  zero,  the  constitutive  equation  for  the  axial  plastic  strain  rate 
from(2.4b)  implies  g  =  0  which  is  also  a  contradiction.  Hence  we  can  only 
have  g  <  0  if  Q  =  0. 
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Figures 


Comparison  of  the  stress  trajectory  in  plane  (measured  in 

MPa),  between  the  theoretical  prediction  ( — )  and  the  exper¬ 
imental  data  (-  -  -)  from  Lamba  and  Sidebottom  (1978a)  for 
OFHC  copper;  the  numerical  calculation  of  the  theoretical 
trajectory  was  effected  with  the  use  of  (3.4),  (3.15),  (3.16)  and  the 
constants  (4.1).  The  prescribed  strain  path  in  the 
plane  is  shown  in  the  inset.  After  initial  elongation  from 
the  undeformed  state  to  point  Q,  the  nearly  elliptical  path 
QRSTQ  is  repeatedly  traversed.  While  the  path  in  the  e^-e^ 
plane  is  cyclic,  the  remaining  nonzero  components  of  strain 
do  not  necessarily  return  to  the  same  values  at  the  end  of 
each  cycle.  Comparison  with  the  experimental  data  is  shown 
for  the  first  few  cycles  only,  since  the  portion  of  the 
trajectory  for  additional  cycles  would  crowd  the  figure. 

Comparison  of  the  theoretical  ( — )  axial  stress  response, 
measured  in  MPa,  with  the  corresponding  experimental  data 
(-  -  -)  of  Lamba  and  Sidebottom  (1978a)  for  OFHC  copper 
according  to  the  prescribed  strain  path  in  the  inset  of  Fig.  1(a). 
Comparison  with  the  experimental  data  is  shown  for  the  first 
few  cycles  only,  since  the  portion  of  the  curve  for  additional 
cycles  would  crowd  the  figure. 


Captions  for  Figures  (continued) 

Fig.  1(c)  Comparison  of  the  theoretical  ( — )  shear  stress  response, 
measured  in  MPa,  with  the  corresponding  experimental  data 
(-  — )  of  Lamba  and  Sidebottom  (1978a)  for  OFHC  copper 
according  to  the  prescribed  strain  path  in  the  inset  of 
Fig.  1(a).  Comparison  with  the  experimental  data  is  shown 
for  the  first  few  cycles  only,  since  the  portion  of  the  curve 
for  additional  cycles  would  crowd  the  figures. 


Fig.  2(a)  The  theoretical  stress  trajectory  in  s^^-s^  plane  (measured  in 

MPa)  calculated  from  (3.4),  (3.15)  and  (3.16)  for  OFHC  copper  using  the 
constants  (4.1).  The  prescribed  strain  path  in  the  e^-e^  p^ane  is 
shown  in  the  inset.  Calculation  for  a  state  of  saturation  was 
obtained  by  first  cycling  in  shear  along  B'D'B'  in  the  e^-e-^ 
plane  for  four  cycles  followed  by  cycling  along  the  elliptical 
path  B'CD’EB'  for  four  cycles.  After  this  saturation  was  attained, 
the  smaller  path  BDBFB  was  repeatedly  traversed.  Only  the  portion 
of  the  stress  trajectory  which  corresponds  to  the  (post  saturation) 
strain  path  BOBFB  is  shown. 

Fig.  2(b)  Comparison  of  the  theoretical  ( — }  axial  stress  response, 

measured  in  MPa,  with  the  corresponding  experimental  data  (-  -  -) 
of  Lamba  and  Sidebottom  (1978a)  for  OFHC  copper  according  to  the 
prescribed  strain  path  in  the  inset  of  Fig.  2(a).  Only  the 
portion  of  the  stress  response  which  corresponds  to  the  (post 
saturation)  strain  path  BDBFB  is  shown. 


Captions  for  Figures  (continued) 

Fig.  2(c)  Comparison  of  the  theoretical  ( — )  shear  stress  response, 

measured  in  MPa,  with  the  corresponding  experimental  data  (-  — ) 
of  Lamba  and  Sidebottom  (1978a)  for  OFHC  copper  according  to 
the  prescribed  strain  path  in  the  inset  of  Fig.  2(a).  Only  the 
portion  of  the  stress  response  which  corresponds  to  the  (post 
saturation)  strain  path  BDBFB  is  shown. 

Fig.  3(a)  Comparison  of  the  stress  trajectory  in  s^-s^  plane  (measured 

in  MPa)  between  the  theoretical  prediction  ( - )  and  the 

experimental  data  ( - )  from  Lamba  and  Sidebottom  (1978b)  for 

OFHC  copper;  the  numerical  calculation  of  the  theoretical 
trajectory  was  effected  with  the  use  of  (3.4),  (3,15),  (3.16)  and  the 
constants  (4.1).  The  prescribed  strain  path  in  the  e^-e^ 
plane  is  shown  in  the  inset.  Calculation  for  a  state  of 
saturation  was  obtained  in  a  similar  manner  to  that  in  Fig.  2(a). 
After  the  saturation  was  attained,  the  complex  path  0-1,  1-2,  ...» 

7-8  was  repeatedly  traversed.  Only  the  portion  of  the  stress 
trajectory  which  corresponds  to  the  (post  saturation)  strain  path 
0-8  is  shown. 


Fig.  3(b)  Comparison  of  the  theoretical  ( — )  axial  stress  response, 

measured  in  MPa,  with  the  corresponding  experimental  data  (-  -  -) 
of  Lamba  and  Sidebottom  (1978b)  for  OFHC  copper  according  to  the 
prescribed  strain  path  in  the  inset  of  Fig.  3(a).  Only  the 
portion  of  the  stress  response  which  corresponds  to  the  (post 
saturation)  strain  path  0-8  is  shown. 


Captions  for  Figures  (continued) 

Fig.  3(c)  Comparison  of  the  theoretical  ( — )  shear  stress  response, 

measured  in  MPa,  with  the  corresponding  experimental  data  ( - ) 

of  Lamba  and  Si  debottom  (1978b)  for  OFHC  copper  according  to  the 
prescribed  strain  path  in  the  inset  of  Fig.  3(a).  Only  the 
portion  of  the  stress  response  which  corresponds  to  the  (post 
saturation)  strain  path  0-8  is  shown. 


Fig.  4  A  geometrical  representation  of  the  function  $  (which  characterizes 

strain-hardening)  as  a  surface  exhibiting  its  dependence  on  the 
plastic  strain  components  e^  and  e^2»  plotted  for  fixed  values 
of  the  total  strains  (in  this  figure  taken  to  be  e-j-^e^O)  f0r 
OFHC  copper  having  the  material  constants  specified  in  (4.1). 

For  different  values  of  e^  and  e12>  the  surface  merely  translates 
parallel  to  the  plane  of  (e^,  e^2)  but  does  not  change  its  shape. 

An  example  of  how  the  value  of  4>  may  be  determined  from  this 
plot  is  also  indicated:  At  a  typical  point  on  the  surface,  the 
values  of  the  plastic  strains  e^ ,  e^2  are  known  (the  coordinate 
curves  on  the  surface  are  drawn  at  plastic  strain  intervals  of 
0.0002)  and  the  corresponding  point  in  the  e^-e^2  plane  is 
located.  Then,  the  value  of  $  is  measured  by  comparing  the 
vertical  distance  between  the  point  on  the  plane  and  the  point 
on  the  surface  with  the  scale  on  the  <t  axis. 
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